Divergent four-point dynamic density correlation function of a glassy colloidal 

suspension: a diagrammatic approach 
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We use a recently derived diagrammatic formulation of the dynamics of interacting Brownian 
particles [G. Szamel, J. Chem. Phys. 127, 084515 (2007)] to study a four-point dynamic density 
correlation function. We re-sum a class of diagrams which separate into two disconnected compo- 
nents upon cutting a single propagator. The resulting formula for the four-point correlation function 
can be expressed in terms of three-point functions closely related to the three-point susceptibility 
introduced by Biroli et al. [Phys. Rev. Lett. 97, 195701 (2006)] and the standard two-point corre- 
lation function. The four-point function has a structure very similar to that proposed by Berthier 
and collaborators [Science 310, 1797 (2005), J. Chem. Phys. 126, 184503 (2007)]. It exhibits a 
small wave vector divergence at the mode-coupling transition. 
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Introduction — It has become clear over the past 
decade that upon approaching the glass transition, the 
liquid's dynamics not only get slower but also become 
increasingly heterogeneous []| . Naturally, one of the very 
first questions that arose after the discovery of dynamic 
heterogeneities was that of their spatial extent. Indeed, 
some of the very first simulational studies of dynamic 
heterogeneities tried to estimate their size and showed 
that it increases upon cooling . The discovery of grow- 
ing dynamic heterogeneities posed a challenging problem 
for the mode-coupling theory [3, 0] of the glass transi- 
tion. First, according to the standard formulation of this 
theory 3], the glass transition is supposed to be a small 
spatial scale phenomenon arising from a self-consistent 
caging of individual particles in their first solvation shells. 
Second, the mode-coupling theory relies upon a factor- 
ization approximation for a four-point dynamic density 
correlation function. Thus, it cannot even be expected 
to describe highly non-trivial spatiotemporal dependence 
of various four-point correlation functions that were in- 
troduced to monitor dynamic heterogeneities. 

The above described problem was first addressed by 
Biroh and Bouchaud (BE) 0. Inspired by an earlier 
study [6] of the so-called schematic mode-coupling equa- 
tions, BB argued, using a field theoretical formulation of 
many-particle dynamics, that the mode-coupling theory 
should be understood as a saddle point approximation 
derived from an action functional. Furthermore, they ar- 
gued that a four-point dynamic density correlation func- 
tion could be calculated by inverting a second functional 
derivative of the same functional and that this proce- 
dure corresponds to a re-summation of ladder diagrams. 
They analyzed the convergence of the diagrammatic se- 
ries and showed that their four-point function diverges 
at the transition point of the mode-coupling theory. In 
a very interesting later contribution Biroli, Bouchaud, 
Miyazaki and Reichman (BBMR) l3| used a more tradi- 
tional, projection operator based [3| version of the mode- 



coupling theory and showed that the matrix which de- 
termines the convergence of BB's four-point function is 
exactly the same as the matrix which describes long 
wavelength properties of a certain three-point dynamic 
susceptibility describing the response of the intermedi- 
ate scattering function to an external potential. Finally, 
in a pair of recent papers [1, 0] Berthier et al. used a 
field theoretical approach to argue that the most diver- 
gent part of the four-point dynamic density correlation 
function is given by, roughly speaking, a product of two 
three-point dynamic susceptibilities. According to Ref. 
8] these three-point susceptibilities can be expressed in 
terms of ladder diagrams (in qualitative agreement with 
BBMR) and thus the most divergent part of the four- 
point function is represented by, r oug hly speaking, the 
sum of "squared ladder" diagrams [lO| . 

The goal of this Letter is to analyze a four-point dy- 
namic density correlation function of a glassy colloidal 
suspension. The reason for considering a colloidal sys- 
tem is twofold. First, experiments provide a wealth of 
information about the motion of individual colloidal par- 
ticles and thus enable detailed tests of various theoreti- 
cal predictions. Second, the simplest model of a colloidal 
suspension, a system of interacting Brownian particles, is 
technically simpler to study than a simple fluid. To ana- 
lyze the four-point function we will use a recently derived 



llj diagrammatic formulation of the dynamics of inter- 



acting Brownian particles. This formulation can use the 
same tools, e.g. re-summations of classes of diagrams, as 
a field theory based approach. It is, however, closer to 
the standard, projection operator based derivation of the 
mode-coupling theory. 

In contrast to BB, we focus on singly connected di- 
agrams, i.e. diagrams which separate into two discon- 
nected components upon cutting through a single propa- 
gator. We will show that a sum of such diagrams results 
in a contribution to the four-point correlation function 
that diverges at the transition point of the mode-coupling 
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theory. The structure of this contribution is very similar 
to the structure of the most divergent part of the four- 
point function obtained by Berthier et al. [1, . 

Since the derivation of the main result is rather tedious, 
we will first define our four-point function and present the 
main result. We will outline the derivation of this result 
in the latter part of this Letter. 

Four-point correlation function — A number of differ- 
ent functions have been introduced to monitor dynamic 
heterogeneities [l3|. In particular, two recent studies 14, 
H, [3 used iV(e-^''-('"iW-'"i)e'i'-('"2W-'"2M(r-ri2)), 
where ri{t) denotes the position of particle i at time t, 
ri{Q) = Ti, ri2 = ri— r2, and N is the number of particles. 
This function quantifies correlations between evolutions 
of two particles that were initially separated by |r|. Its 
Fourier transform with respect to r can be written as 

^g-»k-(ri(t)-r2(t))g«(k-q)-ri2^^ ^^iis work, we cousidcr 

a collective version of the latter expression, 

54(k;q;0 = ^ (n2(k, -k; t)n;(k - q, -k + q)) , (1) 

where n2(ki, k2; t) denotes a projected two-particle den- 
sity [H, 



n2(ki,k2;t) - (l-Po-Vi) 



-iki -Ti (t) — ik2 -Fj (t) 



(2) 

In Eq. ([2]), Vq and Vi are projection operators onto the 
zero-particle and one-particle density, respectively [ll| . 
In the following, we will often use an abbreviated nota- 
tion and, e.g., write n2(ki, k2; i) as ri2(l, 2; t) or as n2{t); 
in addition, n2(t = 0) = n2. 

The projection operators in Eq. Q remove a non- 
zero average value and a non-vanishing projection onto 
the one-particle density. The terms projected out do not 
contribute to the divergent part of 54. Moreover, consis- 
tent usage of projected many-particle densities analogous 
to n2 has important technical advantages 11 , 3| ■ In par- 



ticular, bare inter-particle interactions are automatically 
renormalized by equilibrium correlation functions. 



Contribution to 84 due to singly connected diagrams 
— In this work we focus on the diagrams which separate 
into two disconnected components upon cutting a single 
bond, with each component containing at least a single 
line of bonds from the initial time to the final time t (the 
precise correspondence between S4 and the diagrams will 
be detailed below). The simplest diagram of this type 
is shown on the left in Fig. [TJ In this diagram, bond 
— «- represents a bare propagator and vertices -<, and 
>- represent the left and right vertices, respectively 11|. 
Upon a re-summation, the bare singly connected diagram 
on the left in Fig. [1] is replaced by the diagram on the 
right of this figure. The latter diagram involves two new 
functions, closely related to the three-point susceptibility 
of BBMR, which are represented by T and d , 
and full propagator represented by 



2: 



FIG. 1: Upon re-summation the simplest bare singly con- 
nected diagram on the left is replaced by a renormalized dia- 
gram on the right. 



The main result of this note is the following formula 
for the part of 84 which becomes divergent at the mode- 
coupling transition; this formula results from an approx- 
imate re-summation indicated in Fig. [TJ 



Siik; q;t) 



duj 



nS{q) J 27r^-i 
xG(9;w)xl_q(k;t;w) 



(3) 



Here G{q: uj) is a Fourier transform of the full propagator, 
G{q;uj) = / dte'^'^^G{q;t), and G{k]t) is defined in terms 
of the intermediate scattering function F{k; t), G{k; t) — 
e{t)F{k;t)/S{k). Furthermore, in Eq. © Xq(k; and 
Xq(k; t] Lo) are three-point functions. They are related by 



X^(k;i;-w) 



*Xq(k - q; t; uj). Function Xq(k; t\ uj) 



satisfies the following equation of motion: 



dt' {6{t - t') + M'" (fc; t - t')) '- + ^Xq(k; t-uj)+ dt' I 



dt' Sik)'^''' ' ' ' 7o J (2^)^ |k + q| 

x^;l,(k^k~kOx^(kV-^^a;)f(|k^k^^-0^k+q(k'+q,k-kO ^^^''^^"^,'^''^'^ ^ S^{k;t;uj) (4) 



In Eq. ([5]), Do is the diffusion coefficient of an isolated 
Brownian particle, S{k) denotes the structure factor, n is 
the density and M'^"{k;t) is the irreducible [l3| memory 
function. Moreover, ■i;k(ki,k2) = k • (c(fci)ki -I- c(fc2)k2) 
where c(fc) is the Fourier transform of the direct cor- 



relation function. The source term at the right-hand- 
side of Eq. dH), Sq{k; t; uj) , can be expressed in terms 
of the intermediate scattering function and equilibrium 
correlation functions. This term is regular and finite in 
the q ^ and uj ^ limits, limq^o;Lj^o '5q(k; t; = 



nDoS{0)k'^c{k)F{k;t). 

As explained below, the re-summation leading to Eqs. 
([3]|4]) involves an approximation which is similar to the 
approximation used in the diagrammatic derivation of 
the mode-coupling theory 



Thus, 



to be consis- 
tent we replace the full G, F and M'" in Eqs. dH 
by their approximate forms calculated using the mode- 
coupling equations. Specifically, we assume the standard 
mode-coupling relation which expresses Af"'' in terms 
of F's, M'"{k-t) = (nDo/2)/(dk'/(27r)3)w2(k',k - 
k')F(fc';t)F(lk-k'|;t). 

At this point we note that evolution operator on the 
left-hand-side of Eq. ((4]) with the mode-coupling ex- 
pression for M"^ is exactly the same as the overdamped 
version of the evolution operator derived by BBMR for 
their three-point susceptibility. Indeed, the only differ- 
ence between Eq. (|4]) and the overdamped version of the 
equation of motion (2) of Ref. 01 is the source term. 
Moreover, BBMR emphasized that all of their qualita- 
tive conclusions are independent of the precise structure 
of the source term. Thus, we conjecture that our three- 
point functions Xq(k;t;a;) and Xq(k;ijw) have the same 
scaling behavior as Xq(k;i) of Ref. [7|. In particular, 
these functions exhibit a small wave vector divergence 
upon approaching the mode-coupling transition. 

Furthermore, we note that Eq. ([3]) has the same gen- 
eral structure as the formula derived by Berthier et al. 
(see Eq. (56) of Ref. [1]) Thus, we are lead to 

the same conclusion: the divergence of the four-point 
dynamic density correlation function is intimately con- 
nected to and in fact driven by the divergence of the 
three-point function. 

Derivation — We consider the time-dependent cor- 
relation function of the projected two-particle density, 
(^2(1, 2; t)n2(3, 4)}. For a specific set of wave vec- 
tors this function is proportional to ^4: in the ther- 
modynamic limit we get (ri2(k, — k; i)n2(k — q, k')) = 
n'^Si{\<i- q; t){2TTf5Ck q - k'). 



Following Ref. 



it is possible to derive a hierarchy 
of equations describing the time evolution of (712(^)^2). 
This hierarchy can be rewritten as a hierarchy of integral 
equations which can be solved by iteration. The solution 
can be represented in terms of diagrams consisting of 
the bare propagator, Go{k;t) — 9{t) exTp{—Dok^t/ S{k)), 
represented by — ^, and three vertices, the left vertex, 
V12, represented by -<, the right vertex, V21, represented 
by >-, and the four- leg vertex, V22, represented by X. 
The vertices represent renormalized inter-particle inter- 
actions; they involve the direct correlation function and 
the structure factor rather than a bare interaction poten- 
tial (see Eqs. (55-57) of Ref. 



Ul) 



The diagrams contributing to {n2(t)n2) separate into 
two classes: disconnected and connected ones. For q ^ 
only the latter diagrams contribute to 5'4(k;q;t). Here, 
we only consider singly connected diagrams. A singly 
connected diagram is a connected diagram which sepa- 





FIG. 2: The first few diagrams in the diagrammatic series for 
(712 (ki , k2 ; t)n2 (ks, k4))''™ . 



rates into two disconnected components upon cutting a 
single Go bond, with each component containing at least 
a single line of bonds from the initial time to the final 
time t. The first few such diagrams are shown in Fig. [21 
It is convenient to formulate the diagrammatic expan- 
sion for the singly connected part of {n2{t)n2) divided by 
n'^ S {3) S (A) . For brevity, we will denote this function by 
{n2{t)n2y^"' . In the diagrams contributing to it, we refer 
to the leftmost and the rightmost bare propagators as 
the left roots and the right roots, respectively, and to the 
other bare propagators as bonds. The roots are labeled 
by wave vectors and the bonds are unlabeled. It can be 
showed that {n2{t)n2)^'^'^ is the sum of all topologically 
different singly connected diagrams with two left 

roots, two right roots, Gq bonds, and V12, V21, and V22 



vertices. To evaluate a diagram one assigns wave vectors 
to bonds and integrates over these wave vectors (with a 
{2tt)~^ factor for each integration), integrates over all in- 
termediate times, and divides the result by a symmetry 
number 11] of the diagram; diagrams with odd and even 
numbers of V22 vertices contribute with overall negative 
and positive sings, respectively. 

To re-sum singly connected diagrams we first note that 
the sum of all the "connecting parts" in diagrams in Fig. 
Ogives the full propagator G{q;t) [2l'|. The sums of the 
two remaining parts of singly connected diagrams are re- 
lated to each other by a symmetry operation and thus 
we have to evaluate only one of them. 

We define A"3(ki, q; ka; i, t') as the sum of all con- 
nected, topologically different diagrams with one left 
root, one right root and one left side root (the left side 
root originates from a left "dangling" end of the right 
vertex, V21, or from a left "dangling" end of four-leg ver- 
tex, V22)- The first few diagrams contributing to are 
shown in Fig. [3] 

The last step in the re-summation is the derivation 
of a self-consistent equation for function shown in 
Fig. [4] In this figure, "^^^ represents function X^, a 
thick bond — « — represents the full propagator G(/c; t) 
and represents the reducible memory matrix [ll| . 
Here, we will only argue that equation in Fig. 3] is plau- 
sible [21} . Re-summation of diagrams similar to the first 
and second diagrams in Fig. [3] results in first diagram at 
the right-hand-side of diagrammatic equation in Fig. H) 
Re-summation of diagrams similar to the third diagram 
in Fig. [3]results in the second diagram at the right-hand- 
side of Fig. m Re-summation of diagrams similar to the 
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FIG. 3: The first few diagrams in tfie diagrammatic expansion 
for X3 (ki, q; ka; t'). Note tliat tlie third, fifth and sixth 
diagrams, which contain a single V22 vertex, contribute with 
a negative sign. 




hibits small wave vector divergence at the mode-coupling 
transition. 

Conclusions — We have showed that the contribution 
to the four-point dynamic density correlation function 
of a glassy colloidal suspension due to diagrams which 
separate into two disconnected components upon cut- 
ting a single propagator can be expressed in terms of 
three-point functions which diverge at small wave vec- 
tors. Thus, the four-point correlation function of a glassy 
colloidal suspension exhibits a small wave vector diver- 
gence. In addition, we have derived an explicit formula 
for the dominant part of the four-point function. 

Acknowledgments — I gratefully acknowledge the sup- 
port of NSF Grant No. CHE 0517709. 




FIG. 4: Diagrammatic self-consistent equation for function 
Xi(ki,q; k2;t,t'). 



fourth diagram in Fig. [3] results in the third diagram 
at the right-hand-side of Fig. [H Re-summations of dia- 
grams similar to the fifth and sixth diagrams in Fig. [3] 
result in the fourth and fifth diagrams at the right-hand- 
side of Fig. m respectively. Diagrams similar to the 
seventh diagram in Fig. [3] are neglected. As mentioned 
above, diagrams neglected in our re-summation resem- 
ble diagrams which are neglected in the approximation 
leading to the mode-coupling expression for M" 
Parenthetically, diagrams resulting from the 
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self- 
consistent equation showed in Fig. |4] could be written as 
ladder diagrams. However, these ladder diagrams would 
be quite unusual in that their beams would correspond 
to time running in opposite directions. 

Next, we define three-point function Xq(k;i;i') in 
terms of X^. X^k, q; k'; i, t') = Xq(k; i; i')(27r)^5(k -|- 
q— k')/5(fc'). We introduce its frequency dependent ver- 
sion Xq(k; t; Lv) — /J dt'e*"*'xq(k; t; t'). It can be showed 
that the self-consistent equation for showed in Fig. 
H] is equivalent to the equation of motion ([4]) for the fre- 
quency dependent function Xq(k;t;w). A symmetry re- 
lation mentioned above results in the relation between 
and x^ given below Eq. Finally, we express the 
sum of singly connected diagrams in terms of G, and 
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